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Abstract —The pilot reuse issue in massive multi-input multi¬ 
output (MIMO) antenna systems with interfering cells is closely 
examined. This paper considers scenarios where the ratio of the 
channel coherence time to the number of users in a cell may he 
sufficiently large. One such practical scenario arises when the 
number of users per unit coverage area cannot grow freely while 
user mobility is low, as in indoor networks. Another important 
scenario is when the service provider is interested in maximizing 
the snm rate over a fixed, selected number of users rather than 
the snm rate over all users in the cell. A sum-rate comparison 
analysis shows that in such scenarios less aggressive reuse of pilots 
involving allocation of additional pilots for interfering users yields 
significant performance advantage relative to the case where all 
cells reuse the same pilot set. For a given ratio of the normalized 
coherence time Interval to the number of users per cell, the 
optimal pilot assignment strategy is revealed via a closed-form 
solution and the resulting net sum-rate is compared with that of 
the full pilot reuse. 

I. Introduction 

The MIMO antenna technology has become an essen¬ 
tial component of modern cellular wireless communication 
systems. An exciting recent development related to MIMO, 
collectively referred to as massive MIMO or a large-scale 
antenna array system, is the potential deployment of a very 
large number of antenna elements at base stations (BSs). The 
pioneering work of fH has demonstrated that in time-division 
duplex (TDD) operation with uplink training for attaining 
channel state information (CSI), the effects of uncorrelated 
noise and fast fading disappear as the number of BS antennas 
increases without limit, while only the effect of degraded 
channel estimation due to pilot contamination from the reuse 
of the same pilot set in interfering cells limits the sum-rate 
performance. 

Systematic methods to mitigate pilot contamination effect 
have been studied by various researchers la El El. Shifting of 
pilot frames corresponding to neighboring cells was discussed 
in El, while El suggested coordinated pilot assignment using 
second order statistics for the channels. In El, efficient outer 
multi-cellular precoding under the assumption of appropriate 
cooperation among BSs was introduced for mitigating the pilot 
contamination effect. All known works related to combating 
pilot contamination deal with improving channel estimation 
and rarely mention the potential of allowing more orthogonal 
pilots. 


In Appendix E of IT], Marzetta briefly discusses the per¬ 
formance impact of serving a smaller number of terminals 
than the number of available orthogonal pilot patterns. He 
concludes that while the per-user signal-to-interference ratio 
(SIR) would improve due to reduced pilot contamination (as 
some pilot patterns can be assigned to interferers), the per- 
cell throughput will suffer as the number of users has been 
reduced. This observation is based on the fact that the mean 
throughput per cell improves as the logarithm of SIR while it 
degrades more rapidly as a linear function of the decreasing 
number of users. Let Tcoh be the coherence time interval 
and Tdei be the channel delay spread. It is convenient to 
express the coherence time interval as a dimensionless quantity 
Ncoh = Tcoh/Tdei, after normalization by the channel delay 
spread. Let K denote the number of users in a cell. If the sum 
rate is predominantly a linear function of K, then the net sum 
rate, after discounting the time allocated for the pilot, becomes 
proportional to K{1 — K/Ncoh) assuming K orthogonal pilots 
are used to train K channels. It is then easy to see that the net 
rate is maximized when K = Ncoh/‘I" Marzetta’s message on 
pilot reuse is basically that given a fixed Ncoh, the sum rate will 
be maximized when K is set to Ncoh/"^ ™d the corresponding 
Ncoh/‘^ pilots should be fully reused across all cells; reducing 
K below Ncoh/‘I, in an attempt to reduce interference will only 
hurt the overall throughput. 

This message is certainly correct but is not applicable to 
cases where K cannot simply be adjusted for a given Ncoh- In 
particular, we are concerned with the scenarios where the ratio 
Ncoh/K is considerably larger than 2. Such a scenario arises 
when user mobility is low and at the same time the number of 
users that can be accommodated in a cell cannot be increased 
arbitrarily, as in indoor networks. Another feasible scenario is 
when the service provider wishes to maximize the sum rate 
over a limited number of users (e.g., high-paying customers 
demanding guaranteed data rates) rather than over all users in 
the cell. In fact, higher Ncoh/K ratios due to a relatively small 
number of users per cell reflect meaningful situations in 5G in 
light of the trend toward smaller cells. 

In this paper, we evaluate the maximum net sum rate for 
any given ratio of Ncoh/K and And the corresponding optimal 
pilot reuse strategy under some mild constraints on cell geom¬ 
etry and assuming an infinitely large BS antennas. We show 
that as the ratio Ncoh/K increases, it is beneficial to actually 
prevent full pilot reuse across neighboring cells, i.e., it makes 


sense to allow additional pilot patterns to remove interference 
from users in the neighboring cells even at the expense of 
increased training time. The optimal pilot assignment solution 
is in general not trivial and requires a careful mathematical 
derivation. 

A simple example illustrates when it becomes sensible 
to depart from the full pilot reuse. Using the well-known 
asymptotic result for a large number of BS antennas HI and 
accounting for the pilot portion of the coherent interval not 
used for the actual data transmission, the net user information 
rate for the full pilot reuse (or a reuse factor of one) is given 
by 

Cnet,l = (1 - K/N,oh) log 2 (l + SIR). (1) 

Assume that the interference is coming from a single user 
located at the nearest vertex of a neighboring cell. The cells are 
hexagonal. For comparison, now assume a pilot reuse factor 
of three, where the entire cells are divided into three equal 
partitions, just as in the familiar case of the frequency reuse 
factor three. Again assume for simplicity that there is only one 
interferer who is located in the nearest vertex of the nearest 
cell using the common set of pilots. Given that only the cells 
in a given partition reuse the same pilot set, the distance of the 
worst case interferer from the BS is doubled compared to the 
case for a reuse factor of one. This improves the SIR by the 
factor (2'*^)^ = 4'^, where 7 represents the exponent of signal 
decay as a function of distance. On the other hand, the pilot 
interval is now increased by a factor of three to accommodate 
three distinct orthogonal pilot sets. Accordingly, the net user 
rate for the pilot reuse factor 3 is 

Cnet,3 = (1 - ^K/N,oh) log 2 (l + 4'^SIR). (2) 

It is straightforward to show that Cnet ,3 becomes larger than 
Cnet,i, when 

iV^^3+ log,S/ii 

K 7 

With 7 = 3.8 and an SIR of, say, 3.7 dB, we observe that as 
long as Ncoh/K > 6.2, it is better to go with the pilot reuse 
factor 3 relative to the full reuse. 

The exact threshold value for Ncoh /K for which less than 
full pilot reuse becomes beneficial will change depending on 
the assumed locations of the user and interferers as well as the 
cell partition structure, but it remains true that for given SIR, 
coherence time and number of users in a cell, there exists an 
optimal value of the pilot reuse factor that is not equal to one 
in general. We also show that there exists an optimal way of 
partitioning cells in terms of pilot reuse; we provide a closed- 
form solution to the problem of optimal partitioning wherein 
only the cells in a given partition reuse the same set of pilots. 


II. Pilot Contamination Effect for Massive MIMO 
Multi-Cellular System 


A. System Model 

We assume that the network consists of L hexagonal cells 
with K users per cell who are uniform-randomly located. 
Downlink CSIs are estimated at each BS by uplink pilot 
training assuming channel reciprocity in TDD operation. The 
channel model of HI is assumed in this paper. The complex 
propagation coefficient p of a link can be decomposed into 
a complex fast fading factor h and a slow fading factor 13. 
Therefore, the channel between the BS antenna of the 
cell and the user of the cell are modeled as: 

dmjkl — hmjkl s/I3jkl • (4) 


The slow fading factor, which accounts for the geometric 
attenuation and shadow fading, is modeled as 


f^jkl 


^jkl 


(5) 


where Vjki is the distance between the user in the cell 
and the base station in the cell. The parameter 7 represents 
the signal decay exponent, while Zji-i is a log-normal random 
variable. 


B. Pilot Contamination Effect 


The pilot contamination effect is the most serious issue 
that arises in multi-cell TDD systems with very large BS 
antenna arrays. For uplink training, each BS collects pilot 
sequences sent by its users. Usually, orthogonal pilot sequences 
are assigned to users in a cell so that the channel estimate 
for each user does not suffer from interferences from other 
users in the same cell. However, the use of the same pilot 
sequences for users in other cells cause the channel estimates 
to be contaminated, and this effect, called pilot contamination, 
limits the achievable rate, as M, the number of BS antennas, 
tends to infinity. 

According to HI. under the assumption of a single user per 
cell, the achievable rate during downlink data transmission for 
the user in the cell contaminated by users with the same 
pilot on other cells is given for a large M by 


l0g2 





(6) 


where ffi is the slow fading component of the channel between 
the BS and the interfering user in the cell. In the limit 
of large M, the achievable rate depends only on the ratio of 
the signal to interference due to the pilot reuse. 


This paper is organized as follows. Section II gives an 
overview of the system model for massive MIMO and the pilot 
contamination effect. Section III discusses the assumptions 
made on cell geometry and basic partitioning steps needed 
to utilize longer pilots, and presents the mathematical analysis 
for specifying the optimal pilot assignment strategy for the 
case of single user per cell {K = 1). Simulation results are 
also presented in this section. In Section IV, the solutions 
are extended to the general case with multiple users per cell. 
Section V draws conclusions. 


III. Pilot Assignment Strategy for Single User 
Case 

In this section, we provide analysis on how much time 
should be allocated for channel training given a coherence time 
and how the pilot sequences should be assigned to users on 
multiple cells. We derive optimal pilot assignment strategy, 
which mitigates the pilot contamination effect and maximizes 
total bits transmitted in a given Tcoh- Our analysis considers 
using pilot sequences possibly longer than the number of users 






Fig. 1. 3-way Partitioning 


Fig. 3. Example of coloring, tree structure, and pilot assignment vector 
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Fig. 2. Hierarchical set pailitioning 


in each cell, while orthogonality of the pilots within a cell is 
guaranteed. For simplicity, all cells are restricted to have only 
one user, i.e, K — 1 in this section. The results are extended 
to iT > 1 in the next section. 

A. Hexagonal-Lattice Based Cell Clustering 

Consider L hexagonal cells. Imagine partitioning these 
cells into three equi-distance subsets maintaining the same lat¬ 
tice structure as depicted in Fig.[T] This partitioning is identical 
to the familiar partitioning of contiguous hexagonal cells for 
utilizing three frequency bands according to a frequency reuse 
factor of three. 

It can easily be seen that each coset, having the same 
hexagonal lattice structure, can be further partitioned in the 
similar way. The partitioning can clearly be applied in a 
successive fashion, giving rise to the possibility of hierarchi¬ 
cal set partitioning of the entire cells. In the tree structure 
of Fig|^ the root node represents the original group of L 
contiguous hexagonal cells, and the three children labeled 
L/3 corresponding to the three colored-cosets of Fig. Also, 
applying a 3-way partitioning to a coset results in additional 
three child nodes with labels L/9. Note that a node at depth 
i corresponds to a subset of L3“* cells. A coloring scheme is 
also employed to represent the same pilot sets with common 
color. The number of leaves with different (non-white) colors 
represents the number of different orthogonal pilot sets used. A 
leaf (end node) with a single cell would correspond to a depth 
of log 3 L, but we do not allow such leaves in our analysis as 
this means there would be users with no pilot contamination, 
thus driving the average per-cell throughput of the network 
to infinity as M grows. This particular situation would not 


lend itself to a meaningful mathematical analysis. Thus, the 
maximum depth of a leaf in our tree is set to log 3 L — 1. 

B. Pilot Assignment Vector 

The pilot assignment method can conveniently be formu¬ 
lated in a vector form. Let p be a vector with element pi, 
i — 0, • • • , log 3 L — 1, representing the number of leaves at 
depth i of the partitioning tree. For example, for the tree of 
Fig. we have p = (0,2, 3,0), as there are two leaves at 
depth 1 and three at depth 2. 

Definition; For the given L cells with Ff = 1, is a set of 
valid pilot assignment vectors based on 3-way partitioning and 
is given by 

-Pl = {p= (Po,Pi,'-- ,Fiog 3 L-i) I 0<pi<3* 

logs 

Pi: integer and ^ = 1} (7) 

Definition; For the given L cells, Npii{p) is the length of the 
valid pilot assignment vector p = {po,Pi, - ■ ■ ^Piog^ l-i) and 
is given by 

logs ^-1 

Npilip) = ^ Pi (8) 

i=0 

The pilot length Npii{p) represents the number of orthogonal 
pilots corresponding to the given pilot assignment vector p. As 
an example, for the pilot assignment strategy shown in Fig. 
we have Npii{p) = 5. 

Notice that the users in different cells belonging to a 
given leaf experience pilot contamination. The severity of the 
contamination depends on the depth of the leaf. Every time the 
depth is increased, the distance between interfering cells (cells 
that reuse the pilot set) increases by a certain factor, as can 
be observed from Fig. In fact, the distance grows geometri¬ 
cally as the depth increases. According to (j^, the achievable 
rate is determined by the (3 values of the interfering users, 
which in turn depend on the distances of the interferers from 
the BS. The distance growth manifests itself as the reduced 
pilot contamination effect or an improved SIR, increasing the 
throughput. More specifically, the throughput grows roughly 
with log 2 (Q!^'i'S'/i?i) = 27 log 2 (Q;) -|-log 2 SIRi, where SIRi 
is the reference SIR corresponding to the pilot reuse factor 1 


and a is the parameter that represents the distance growth. It 
is clear that a geometric growth of a gives rise to a linear 
increase in the throughput. Letting Ci be the rate of a user at 
depth i, this is to say that Ci increase linearly with depth i. 

The per-cell sum rate for the network can be expressed as 

^ loggL-l logjL-l 

CsumiP) = J^ Y. Li-^P^C^= Y 3-VzC,. (9) 

2—0 2 — 0 

The per-cell net sum rate, accounting for the fact that useful 
data gets transmitted only over the portion of the coherence 
time not allocated to the pilots, is given by 

C^etiP.N^ou) = (10) 


transition element to = 1 triggers a (3-way) partitioning at 
depth 0, temporarily creating a pilot vector (0,3,0,0). The next 
element = 1 induces a (3-way) partitioning on one of the 
3 existing partitions at depth 1, thereby giving rise to a new 
pilot vector (0,2,3,0). Since the next transition vector element 
is zero, the partitioning stops. The transition elements also 
point to the number of white nodes at each depth, as can be 
confirmed in Fig. 

Lemma 2. Any transition vector t = i^iog 3 L- 2 ) 

originated from p S n{NpQ) satisfies 

fO<L<3* 0<i<log3L-2 

Npo-1 (16) 


logs L-2 

Y = 


We shall use Cnet as the objective function for finding optimal 
pilot assignment strategies. 

C. Closed-Form Solution of Optimal Pilot Assignment Strategy 

The optimal pilot assignment vector Popt{Ncoh) for the 
given normalized coherence time Ncoh and the number of cells 
L is: 

Popt{Ncoh) = argmax Cnet{p,Ncoh)- (H) 

P&Pl 

We also note that for K = I, all valid pilot assignments yield 
odd pilot lengths, as formally stated in Lemma [T] 

Lemma 1. For K = 1 and given L, {Npifip) : p S Pl} = 

,f}. 

Before giving the first main theorem, it is useful to define 
the pilot assignment vector that maximizes the per-cell sum 
rate Csum with a finite pilot length constraint: 

Popt(^po) = argmax C sum (P) (12) 

pGn(A^po) 

where 

^{Npo) = {p€PL\Npa{p) = Npo}. ( 13 ) 

We will first find a closed-form solution for p^^j and then find 
eventually Popt by exploring its relationship with the former. 

Definition: For each valid pilot assignment vector p = 
(pO)Pi) • • • jPiogsL-i), the corresponding transition vector 
t = (fo, fi, • • • , flogs t^- 2 ) is defined as: 

r to = 1 - po 

\ti = -pi -\- iti-l. 1 < J < logg L-2 
The inverse relationship exists: 

Po = 1 - fo 

Pi = 3ti_i — ti 1 < * < logo L — 2 

PlogsL-l = 3flogsL-2- 

The first two equations of ( |T5] l come from ( [T4) i, while the 
last equation is from the fact p,Yi^=Y~^ — 1- Each 

element of t represents the number of transitions at depth 
i from the conventional (i.e., full pilot reuse) assignment 
(1,0,-•• ,0). For example, p = (1,0,0,0) turns to p = 
(0,2,3,0) via a transition vector t = (1,1,0). The first 


(14) 


(15) 


We further define the index function xi^po) that identifies 
the first non-zero position of Popt(lVpo): 

^ AT _ 1 

x{Npo) = min{fc | X! i 

2^0 

Theorem 1, For K = 1, given L and NpQ, the optimal 
pilot assignment vector p'^p^{Npo) = (p(,,- • • with 

respect to Csum, has its components as follows: 



Y_r- 


s=0 



0 


NpQ-1 

2 

* = Xi^po) 


i = xiNpo) + 1 

s=0 / 

otherwise 


(18) 


Flere we only provide a brief sketch of the proof due to 
space limitation. The sum rate in (j^ can be expressed using 
the transition vector t via 

logg L—2 

CsuTn{t) = Co+ Y f^3-‘(C,+l - Q). (19) 

z=0 

Because Ci is a linear function of i, the difference (Ci+i — Ci) 
is a constant. Therefore, all we have to do is to find the optimal 
t which maximizes ^ fi3“*. 

Under the condition of Lemma for maximizing 
^■°fg L3“®, we should give the largest values U = 3* 

to the lower indices i < xi^po) ^nd the remaining val¬ 
ues to the rest of the indices. The final solution for t' = 

(foi fli ■ ■ ■ I flogs i- 2 ) 

x{Nj,o) — l 

f, ^ tVgpi _ g 3* 

2^0 

0 

V 

which leads to ( fTS] ) via 

For a given NpQ, more than one valid pilot assignments may 
exist. For example, if L = 81 and NpQ = 7, p = (0,1,6,0) 
and p = (0,2, 2,3) are valid vectors, but Theorem H reveals 
that PoptiNpo) = (0,1,6,0). 


f < Xi^po) 

^ = xiNpo) (20) 

f > xiNpo), 












Optimal pilot assignment (L = 81, A" = 1) 


From Theorem [T] a certain trend relating the optimal 
pilot assignment vectors p'^p^{Npo) and p'^p^{Npo + 2 ) can 
be observed as follows. 

Corollary 1. For two pilot lengths NpQ and NpQ + 2, the two 
corresponding optimal pilot assignment vectors p'^p^^Npo) = 

(p 5. • ■ ■ , Kog 3 L- 1) and p'^pt {Npo + 2) = {p*o *, • • • , p**^^ ^) 

exhibit the following relationship: 

(p* -1 i = x{Npo) 

pr = <p:+3 1 = xiNpo) +1 ( 21 ) 

[p* otherwise 

For a given NpQ, the optimal assignment vectors p'opti^po) 
and p'^p^{NpQ + 2 ) show a predictable pattern of tossing 1 
from the left-most non-zero component to increase the adjacent 
component by 3. For example, in the case of L = 81, the 
optimal assignment Pop((7) = ( 0 , 1 , 6 , 0 ) can be transformed 
by reducing the second component by 1 and increasing the 
third one by 3, which results in the next optimal assignment 
for Npo = 9: Popt(9) = (0,0,9,0). It can be seen that there is 
a tendency to reduce the left-most non-zero values which give 
the most severe pilot contamination. 

We now set out to find Popt{Ncoh)- Theorem and 
Corollary already identify, given the fairly mild constraints 
of hexagonal cells and equi-distance partitioning, the optimal 
pilot assignment strategy maximizing the sum rate for the 
chosen pilot sequence length. The next step is to find the 
relationship between the normalized coherence time Ncoh and 
the optimal pilot sequence length. 

First, write the net sum-rate as 

C'netiPopti^po) 1 ^coh) 

= (22) 

which is an increasing function of Ncoh and crosses the 
horizontal axis once at Ncoh = ^po- Moreover this function 
saturates to Csum{p'opti^po)) for very large Ncoh- Imag¬ 
ine plotting this function for all possible values of NpQ = 
1, 3, 5, 7, • • •. As NpQ increases, the zero-crossing is naturally 
shifted to the right while the saturation value moves up. More 
specifically, the Cnet curve for Npo = 2 n — 1 and that 
for Ncoh = 2n + 1 intersect once. On the left side of this 
intersection point, the Cnet curve for Npo = 2 n — 1 is above 
the latter while on the right side, the latter curve is higher 
than the former. Let the horizontal value of this intersection 
point be Ncoh = A„. Applying the same argument, the two 
Cnet curves for Npo = 2n + 1 and Npo = 2n + 3 intersect 
at A„+i somewhere to the right of the previous intersection 
point An- In between A„ and A„_|_i, the Cnet curve for 
Ncoh = 2n+l yields the highest values, indicating that when 
Ncoh falls between the two intersection points, the optimal 
pilot length is 2n + 1. It can be shown that the intersection 
points are given by 

X(2n-1)-1 \ 

2n-l- -fl, 1 < n < Nl 

i^O I 

(23) 

where ^(n) = ^)C^i^ 2 n-l)/{C^i 2 n-l) + l - C'x( 2 n- 1 )), 

with Ci already defined in Section III-B, and N^ is the 


TABLE I. 


Nooh 

Popf (-^ coh ') 

^ pil { Popt { Ncoh )) 

0 ~ 6 

(1.0.0,0) 

1 

7 ~ 23 

(0.3.0,0) 

3 

24 ~ 27 

(0.2.3,0) 

5 

28 ~ 31 

(0.1.6,0) 

7 

32 ~ 70 

(0.0.9,0) 

9 

71 ~ 74 

(0.0. 8, 3) 

11 

75 ~ 78 

(0.0.7, 6) 

13 

79 ~ 82 

(0.0. 6, 9) 

15 

83 ~ 86 

(0,0,5.12) 

17 

87 ~ 90 

(0,0,4.15) 

19 

91 ~ 94 

(0,0,3,18) 

21 

95 ~ 98 

(0,0,2.21) 

23 

99 ~ 102 

(0,0,1,24) 

25 

103 ~ 

(0,0,0,27) 

27 


number of all possible pilot lengths minus one. Since Npu = 
1, 3, • • • , L/3 from Lemma B we have Nl = - !)■ 

We now state our second main theorem without a formal 
proof (which will be given elsewhere). 

Theorem 2. For K = 1, given L and Ncoh, if the normalized 
coherent time Ncoh A in between two adjacent time points 
An and A„_|_i, i.e.. An < Ncoh < A„+i, then the optimal 
assignment vector Popt{Ncoh) satisfies 

Npil {Popt{Ncoh'}^ — 2 tI 1 

Popt{Ncoh) = Popt(2tt + 1), (24) 

where p'cp^{2n + 1 ) is the vector maximizing Caum given the 
fixed pilot length (2n -f 1). 

D. Simulation Results 

Simulation is necessary to compute the Ci values in (j^. To 
compute Ci, the terms in ^ need to be generated pseudo- 
randomly according to the assumed underlying statistical prop¬ 
erties. For this, we assume; a signal decay exponent of 7 = 3.8, 
the shadow-fading standard deviation of (ishadow = 8.0 dB, 
and a cell radius of rc = 1600 meters with a cell-hole radius 
of Th — 100 meters, to be consistent with the parameter 
values used in H]. The user locations are uniform-random. 
To generate each Ci value, an average is taken over 100,000 
pseudo-random trials. The number of cells are fixed to L = 81 
and the number of users in each cell in this section is AT = 1 . 

Once the Ci values are computed, the optimal pilot length 
and the pilot assignment vector as well as the net throughput 
can be obtained for various coherence intervals Ncoh - Table |I] 
shows the optimal pilot assignment results for various values 
of Ncoh- We confirm that the simulation results of Table ||] are 
consistent with the mathematical analysis given in Theorem]^ 

Fig. 0 shows the average achievable net rates for various 
pilot assignments versus normalized coherence interval. The 
random assignment means that a pilot sequence is chosen 
randomly and independently from Npifipopt (Ncoh)) orthog¬ 
onal pilots, and assigned to each user. Therefore, the optimal 
assignment and the random assignment use the same amount 
of pilots for any given Ncoh- It can be seen that a substantial 
performance gain is obtained using the optimal method com¬ 
pared to the full pilot reuse case as Ncoh increases beyond 7. 
The random assignment is worse than the full reuse initially 
but eventually outperforms the latter as Ncoh grows. 








Theorem 3. For given NpQ, L and K, the optimal pilot assign¬ 
ment vector p'^p^{Npo,K) = {po,--- ,p[og^L-i) maximizing 
Csum has its elements p'^ in the form: 



TK3*-^ 


t=o 



-K 

2 

* = x{Npo,K) 

i-2 \ 

y:A'3M 

t = x{Npo,K) + l 

/ 

otherwise 

(30) 


We have another generalized definition; 
Definition; 


Fig. 4. Net rate for various pilot assignments (L = 81, K = 1) 


For Ncoh = 20,40 and 80, for example, the optimal 
assignment method has 35%, 59% and 81% higher net rates 
Cnet than the full pilot reuse assignment, respectively. As co¬ 
herence time increases, the benefit of allocating more time for 
pilots is considerable. Also, the non-shrinking performance gap 
between the optimal assignment and the random assignment 
indicate that structured optimal assignment is required for a 
given pilot time, in order to maximize the net throughput of 
the network. 


aw = 


r}{n,K) — l 

2 I 2n - 1 - K3" + K^{n, K) \ + K 

i=0 


K 


for 1 < n < Nk,l, where 


'7j{n,K) = x{2n + K-2,K) 

^{n,K) = - 

LK _ 

— 


At the initial point, = 0. 


(31) 


(32) 


IV. Multi-User Analysis 

We now consider multiple users per cell. Due to space 
limitations we will state our main results in the form of two 
main theorems without proof. Interestingly, the optimal pilot 
assignment vectors take a similar form to those for the K = 1 
case. 

A. Analysis for Multi-User Case 

The set of valid pilot assignment vectors generalized to 
AT > 1 is defined as; 


Using pT] ), Theorem]^ can be generalized as follows; 

Theorem 4. For given L,K and Ncoh, if the ratio Ncoh/K 
is in between two adjacent time points and i.e., 

a1^^ < ^ then the optimal assignment vector 

PoptiNcoh,K) maximizing CnetiPiN^oh) satisfies 

Npii (PoptiNcoh, K))=2n + K 
Popt{Ncoh,K) = p'^p^{2n +K,K), (33) 

where p'^p^(2n -f K,K) is the vector for fixed pilot length 
(2n -\- K) that maximizes Cgum for given L and K. 


Pl,k = {p={po,Pi,--- .Pion^L-i) I 0<p, <a:3* 

logs tj-l 

Pi ', integer and ^ Pi3“* = K} (25) 

i=0 

Also, generalize the following definitions; 

Popt{Ncoh,K) = arg max Cnet{p, Ncoh) (26) 

PGPl.k 

U(Vpo, K) = {p€ Pl,k I Npuip) = Npo} (27) 
p'opt{Npo,K)= arg max C 

sum (p) (28) 

pGO{Npo,K) 

k 

xiNpo, K) = mm{A; | ^ iF3* > (29) 

i=0 

We present the following theorem for general K: 


B. Simulation Result 

In Table [n| the optimal assignment vectors and pilot lengths 
for various Ncoh are shown, assuming L = 81 and K = 2. 
Like in the case for K = 1, the optimal pilot assignment 
vectors have a predictable form. As Ncoh increases, the opti¬ 
mal pilot sequence gradually becomes longer. Moreover, the 
optimal assignment vectors show a pattern of tossing 1 from 
the left most non-zero component to increase the adjacent 
component by 3, for an example (0,6,0,0) —>■ (0,5,3,0) 
(consistent with Corollary [^. 

Fig. shows the plots of the per-user net rates Cnet /K 
versus TVoft, / K. These results are actually found for K = 14, 
but very similar numerical results are obtained for the optimal 
scheme irrespective of the particular values of K while the 
results for the full pilot reuse are identical across all values 
of K. As a case in point, it can be seen that the plots in Fig. 

obtained for Ff = 1 are nearly identical to those in Fig. 
corresponding to K = 14. 
















TABLE II. Optimal pilot assignment information 
(L = S1,K = 2) 


Ncoh 

coh 1 I^) 

Npii{popt{Ncoh, K)) 

3 ~ 11 

(2, 0.0.0) 

2 

12 ~ 15 

(1.3. 0.0) 

4 

16 ~ 45 

(0.6. 0.0) 

6 

46 ~ 49 

(0.5. 3.0) 

8 

50 ~ 53 

(0.4. 6.0) 

10 

208 ~ 

(0.0.0.54) 

54 



Fig. 5. Per-user net rate versus 



changes as K decreases while Ncoh is held fixed. Notice that 
the maximum net sum-rate occurs at Ncoh/K = 2 and at 
this point the optimal scheme reduces to full reuse, consistent 
with the Marzetta’s analysis HI . However, the message here is, 
again, that when we do not have a control over Ncoh/K = 2, 
the optimal pilot assignment strategy may have a substantial 
net sum-rate advantage over full pilot reuse. 

To appreciate how large Ncoh/K can be in some real- 
world scenarios, take an indoor office wireless channel with 
Tcoh = 50 microsec and Tdei = 50 nsec, yielding Ncoh = 
1000. If the user density cannot be allowed to be more than 
K = 20 users per cell, we would be focusing on Ncoh/K = 50 
and at this point, the optimal scheme gives a 67% net sum- 
rate improvement over full pilot reuse. As another example, 
consider an urban outdoor environment with a fairly high user 
mobility giving rise to a 1 msec coherence time interval. With 
a 2 microsec delay spread, for example, this gives Ncoh = 
500, and assuming not more than K = 25 users are to be 
served, we are interested in the net sum-rates at Ncoh/K — 20. 
From either Fig. or Fig. we see that a net throughput 
improvement of 35% is possible via optimal pilot assignment, 
relative to full pilot reuse. 

V. Conclusion 

In massive MIMO systems with interfering cells, when 
the coherence time interval Ncoh and the number of users 
K are given, finding the appropriate portion of the coherence 
interval to be allocated to channel training is not trivial. This 
paper has provided an analytical solution to finding the optimal 
training time for any given ratio of Ncoh/K and along the 
way showed that the net throughput could be improved in 
general by allowing neighboring cells to use different sets of 
pilot sequences. Assuming hexagonal cells and equi-distance 
hierarchical partitioning, an optimal pilot assignment strategy 
has been identified that gives substantial throughput advantages 
relative to random pilot assignment or full pilot reuse when the 
ratio Ncoh/K is sufficiently large. Finally, we add that it would 
be interesting to further explore pilot assignment strategies 
when the objective is not about maximizing the sum rate but 
rather on guaranteeing some minimal performance level to 
all users or maximizing a weighted sum rate to prioritize the 
services. 


Fig. 6. CnetlNcoh versus Ncoh/K for different pilot assignment schemes 

In fact, it can be shown that for all valid pilot vectors 
with only one non-zero component, p = {K, 0, • • • , 0), 
(0,376,0,-•• ,0), (0,0,9i6,0,--- ,0), •••, (0,---,0,^), 
the plot of Cnet/K verses Ncoh/K does not change with K. 
For the remaining pilot assignment vectors, the plot changes 
very little across different values of K. As can be seen, using 
the optimal pilot assignment scheme the per user net rate 
improves substantially with increasing Ncoh /K, relative to the 
full reuse scheme. 

Fig. 1^ shows Cnet/Ncoh versus Ncoh/K for optimal, full 
reuse and random pilot assignment schemes. Again, AT = 14 
is used to generate these plots (L = 81 throughout this paper), 
but the plots do not change noticeably for different values 
of K. These plots give an insight into how the net sum rate 


References 

[1] T. L. Marzetta, ’’Noncooperative cellular wireless with unlimited numbers 
of base station antennas,” IEEE Trans. Wireless Commun, vol. 9, no. 11, 
pp. 3590-3600, Nov. 2010. 

[2] K. Appaiah, A. Ashikhmin, and T. L. Marzetta, ’’Pilot contamination 
reduction in multi-user TDD systems,” in Proc. IEEE Int. Conf. Commun. 
(ICC), May 2010, pp. 1-5. 

[3] H. Yin, D. Gesbert, M. C. Filippou, and Y. Liu, ’’Decontaminating pilots 
in massive MIMO systems”, in Proc. IEEE Int. Conf. Commun. (ICC), 
Jun. 2013, pp. 3170-3175. 

[4] A. Ashikhmin and T. L. Mai'zetta, ’’Pilot contamination precoding in 
multi-cell large scale antenna systems,” in Proc. Int. Symp. Inf. Theory 
(ISIT), July 2012, pp. 1137-1141. 




















